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In an N-level quantum mechanical system, the problem of unitary feedback stabilization of mixed 
density operators to periodic orbits admits a natural Lyapunov-based time- varying feedback design. 
A global description of the domain of attraction of the closed-loop system can be provided based on 
a "root-space" -like structure of the space of density operators. This convex set foliates as a complex 
flag manifold where each leaf is identified with the coadjoint orbit of the eigenvalues of the density 
operator. The converging conditions are time-independent but depend from the topology of the 
flag manifold: it is shown that the closed loop must have a number of equilibria at least equal to 
the Euler characteristic of the manifold, thus imposing obstructions of topological nature to global 
£h I stabilizability. 
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INTRODUCTION 



The use of feedback in quantum mechanics is limited by the phenomenon of wavefunction collapse following 
a measurement. In this work the problem is bypassed by considering density operators of quantum ensembles 
■ and completely noninvasive measurements. This allows also to relax the requirement of commutativity of 

\Q ' the measured observables and in fact we shall assume to have a complete knowledge of the density operator 

for all times. Although physically this set up is realistic only for some applications (typically nuclear spin 
ensembles 0H3), it is of widespread use for the purposes of model-based quantum control (often under 
the name "tracking control" Il2l l38j l. as it allows to generate control fields in spite of the high complexity 
of open loop control |lCt ITfjl l34|. Furthermore, while the formulation comes from quantum control, our 
Oh! motivations for this work are mostly mathematical, namely feedback design and convergence analysis for a 

class of bilinear control systems living on a particular family of compact manifolds and evolving isospectrally. 
As the system has a drift term which cannot be canceled without incurring in singularities of the control 
law, the most natural problem formulation is to seek for a stabilizer to the periodic orbit drawn by the drift. 
Rather than studying this problem like an orbital stabilization problem we reformulate and solve it as 
a state tracking problem, thus avoiding the obstruction to semiglobal convergence of a periodic orbit, see 
[3?|. Corollary 1.6 (where it is called stability in the large). In fact, with our feedback design the state will 
converge to the orbit, but the entire orbit is not an invariant set, only a point moving along it is invariant. 
As a matter of fact, by passing to a suitable rotating frame, our time-dependent trajectory tracking problem 
can be reformulated completely in terms of time- varying feedback law for the fixed point of a nonautonomous 
system. The Lyapunov design is essentially of the Jurdjevic-Quinn type |24|. for which the usual LaSalle 
invariance principle is applicable in spite of the time-dependence of the closed loo p, a nd does not differ much 
from what has already been proposed in the literature for wavefunctions 0, l2lL l28t l3f| . 

What is nontrival is to ascertain the convergence of certain initial conditions and to provide a global 
description of the region of attraction. In fact, the sufficient condition used in [24[ to prove asymptotic 
convergence and based on the so-called ad-condition or Jurdjevic-Quinn condition 0, is never verified globally 
for N > 2. This is due to the presence of an abelian subalgebra (Cartan subalgebra) that can never be fully 
spanned by ad-commutators alone. It will be shown, however, that the undesired critical points are not only 
unstable but also repulsive, meaning that the Jurdjevic-Quinn condition or, equivalently, the controllability 
of the linearization along the desired orbit , guarantees convergence for all initial conditions outside the set 
of equilibrium points. To attain a complete and time-independent description of the critical set and thus of 
the domain of attraction, a thorough geometric and topologic characterization of the state manifold of mixed 
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density operators is required. A unitary evolution like that appearing in a Liouville equation is isospectral, 
as the eigenvalues of the density operator form a complete set of invariants. Unlike for a wavefunction, the 
state space has dimension and structure which depend on the multiplicities of such eigenvalues. Since these 
form a flag in dimension N, all complex flag manifolds obtained as homogeneous spaces of U(N) (or SU(N)) 
by the Cartesian products of subgroups of dimensions given by the multiplicities are admissible state spaces 
H EUL Hil- Since the set of density operators up to the imaginary unit "overlaps" with the Lie algebra 
u(N) (or su(N), excluding the constant trace), these complex flag manifolds can also be intended as the 
orbits of the (co)adjoint action of U(N) (or SU(N)) on its Lie algebra. This is a well-studied action and the 
structure of its orbits is well-known |19l l27| : for example a fundamental topological invariant like the Euler 
characteristic acquires the meaning of number of nontrivial possible permutations of the eigenvalues of the 
density operator, see also |8l ll.ll Il7l l39|. For the purposes of stabilizability, this is an important feature, 
because it will be shown that each complex flag manifold has a number of "antipodal" points equal to the 
Euler characteristic, and that these points must be equilibria of the closed loop system. In order to give a 
complete description of the region of attraction, we use the resemblance between the set of density operators 
and the Lie algebra su(N), and a few tools deriving from the root space decomposition of a compact Lie 
algebra, namely its orthogonal decomposition into Cartan subalgebra plus root spaces and the invariance 
properties of the root spaces under certain commutators (like the ad-commutators) Q. This "graph-like" 
approach yields simple, time-independent characterizations of all converging initial conditions for a given 
reference orbit and Hamiltonian. Also the Kalman controllability of the linearization admits an intrinsic 
formulation in these terms. The formalism used gives insight into the problem of choosing reference orbits 
having a large domain of attraction. 

It is known 9, 26], that compact manifolds do not admit a global asymptotically stable equilibrium because 
they are not contractible. This is a topological property and corresponds to a set being homotopy equivalent 
to a poi nt [221 . The region of convergence of an asymptotically stable attractor must be in such a homotopy 
class 0, l.37j. For our complex flag manifolds, it will be shown that the antipodal points represent topological 
obstructions to global stabilizability. 

In order to simplify the treatment, an equivalent real representation of density operators is used throughout, 
given by the so-called coherence vector and correspondingto the vector of expectation values with respect 
to a complete orthonormal set of Hermitian matrices l35l ]. It provides a linear representation of the 

adjoint action occurring in a Liouville equation and it allows to formulate the control system in terms of 
standard bilinear systems on smooth manifolds which are real representations of the complex flag manifolds. 

II. DRIVEN LIOUVILLE- VON NEUMANN EQUATION 

With a given Hamiltonian H = Ha + uHb, —%Ha, —iHb € su(N), u G C°°(R) a control field, one can 
form a Schrodinger equation for the wavefunction (in atomic units, K = 1) 

\^) = -i{H A + uH B )\V), m))eS 2N -\ (1) 

or a Liouville-von Neumann equation for the density operator p 

p = -i[H A + uH B , p], p = pt>0, tr(p) = l, tr(p 2 )<l. (2) 

Eq. J2J holds for a quantum ensemble, hence it is more general than (T|: the two being equivalent only when 
tr(p 2 ) = 1 (i.e., p is a rank-one operator: p = |^)(-0|). 

A. Gell-Mann basis and adjoint representation 

The left hand side of (0) contains a conjugation action on state matrices. In order to deal with it, we 
can use one of the features of the adjoint representation, namely its providing a linear representation of 
one-parameter groups of automorphisms of su( N) (see Appendix EJ , and reformulate (0 as a standard 
bilinear control systems on a suitable manifold. To do that, we use the coherence vector representation of 
p, whose key property is that both H and p are expressed in terms of the same complete orthonormal set 
of Hermitian matrices. Let Aq = — }= 1 n and call A the n-dimensional vector of N x N Gell-Mann matrices 
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[20| . n — N 2 — 1. Since 5u(N) contains traceless skew-hermitian matrices, span{— iX} — su(N). Denote 
with t) the Cartan subalgebra of su(N), i.e., the maximally abelian subalgebra in su(N), dim(h) = N — 1. In 
the Gell-Mann basis A, f) corresponds to the N — 1 diagonal matrices. Denote 6 the vector space such that 
su(N) = f) © 6, with f) _L f in a standard biinvariant su(iV) metric (e.g. the Killing metric). While f) is an 
abelian subalgebra, t is only a vector space. In correspondence of this direct sum, we have the decomposition 
of A into Af, and At so that f) = span{— i\t>} and 6 = span{— iAt}. If Ej£ is the elementary N x N matrix 
having 1 in the (j£) slot and elsewhere, then the matrices A are given by 

{A,,,,, 1 < j < N - 1} = {(£n + • ■ • + Stf - jE j+lij+1 )/y/j(j + 1), l<j<JV-l} (3) 
for the diagonal part, and 

{M,Rjt, Uj</a} = {fe+£(j)/^, l^j<£<7v} (4) 
{Ae,3j/, 1 < i < € < N} = [H-Eji + E tj )/y/2 t 1 < j < £ < ivj (5) 

for the off-diagonal part. Calling = span {— iA^sr j^, — iXt.Q.je}, then we have the further splitting of 6 
into "root spaces" 

«= «i* (6) 

with the following commutation relations (see for instance |3j for the details): 

h tit] = Ijt, (7) 



%h tpq] — 



tjq 



if £ p and j =/= q 
if£ = p 
if j = q 

if j = p and £ = q. 



(8) 



Assume 7?^ is diagonal 



"JV 



+ £ n = 0, 



and nondegenerate, i.e., £j ^ f^, j ^ ^, where the £j are supposed ordered: £1 < £2 < ■ ■ ■ < £n- The £j 6 K 
are the energy levels of the (unforced) system QJ, i.e., the eigenvalues of the stationary Schrodinger equation 
Hji\ipj) = £j\ipj) of eigenstates \ipj) = ej (ej is the elementary basis vector), j — 1, ...,iV. Assume also 
that the transition frequencies are nondegenerate, i.e., that the levels are not equispaced £j — £i 7^ £ p — £ q , 
7^ ipi) J / P / 9- Further standard assumptions are that Hb is off-diagonal and Graph(i7e) 
connected. A stronger assumption we shall need is that Hb enables all transitions among adjacent energy 
levels: tr (Hsij.j+i) 7^ V j = 1, ...,N— 1. Beside connectivity of Graph(ffe), it guarantees that all 
"fundamental root spaces" (see 0) are excited by the dynamics. Since H a is diagonal and traceless, —iHa G 
su(N). A nondegenerate element of the Cartan subalgebra (like Ha above) is called regular. It is called 
strongly regular if in addition it has all nondegenerate transitions. 

In terms of A, Ha = fiA ■ X and Hb — hs ■ A, with h,A, hs G K."- Likewise p — poAo + Q • A, where 
g = — L and the coherence vector g is composed of expectation values along the basis elements JSJ-JS}: 
Qf)d = tr (/?Af,j) G R, j = 1, . . . ,N - 1, Qijstji = tr(pX mje ) G K, = tr(pA tj cj^) 6 I, U j < JV. 

Any density p can be split as p = QqXo + Ptj + Pt, or m correspondence of p = goAo + ■ X^ + Qi ■ 

A { = qo\q + J2i^ 3 <N 0t)j^,j + J2i^ 3 <e^N (Qt,Xjt*t,R,jt + Qt.^.jeXt.^.ji)- The nonzero components of the 
coherence vector q uniquely identify a subset of the tjg,. Denote ft(p) the "support" of p in t, i.e., the set of 
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root spaces "touched" by p: f t (p) = p n 6. Also let J" t (p) = {(j'f) s.t. tr (pt^) ^ 0, 1 < j < £ < N} be the 
corresponding set of index pairs. When € Flip), then (gt t st,je, Qt,3,jl) ^ (0,0). Likewise f (>o) = p H f) 
and Ftjip) = {(j) s.t. g^.j 7^ 0, 1 ^ j < N}. In the following we shall use both symbols p and g for densities 
and we shall refer to "diagonal" and "off-diagonal" g with an obvious abuse of notation. The Hilbert- Schmidt 
norm on density operators induces for g the standard Euclidean norm || • ||: tr (p 2 ) = gl + {{g,g)) = £?o + ||f?|| 2 , 
where we indicate with (( • , • )) the M"-Euclidean inner product. Due to the trace-class constraint, the notion 
of distance between the densities pi and p2 having the same purity is d(p%, p-i) = tr (p 2 ) ~ tr {P1P2), see e.g. 
[3^|. or in terms of g: 



d(gi, Q2) 



el 



Qi 



el 



(01,02)) = \\gi\\ 2 - gUi e[0, tr( P 2 )] 



(9) 



Thanks to the use of the same basis for p and the Hamiltonian, up to the imaginary unit the trajectories 
of |(2Jl can be identified with the adjoint orbits of SU(N) on its Lie algebra, see Section ITTT1 for a thorough 
description. Following Appendix EI we can replace the matrix ODE J2J) with the linear vector ODE 



g={A + uB)g, 



(10) 



where A = — iad# 4 = —ih,A ■ &d\, B = —iadH B = —ihs ■ ad.\ and A,B^z ad 5U (jv) C so(n). The Lie algebra 
ad su (jv) is the adjoint representation of Su(N), hence dim(ad su (^v)) = dim(su(iV)) = n. Therefore, for N > 2, 
ad su (Ar) C so(n). 

The isomorphism su(N) ~ ad su (jv) induces an orthogonal splitting also in the adjoint representation: 



ad 



su(JV) 



adf,©ad{. Obviously relations similar to still hold and A = —ih,A ^-adx,,, B 



For C £ su(N), in components C = J2i^ 3 <n c t)d\,3+J2i^<e^N (ct.XjtM&jt + Ct&jtAi&jt) 
indicate with f(,(C), fe(C) the support of C in, respectively, f), 6, of indices !Ft,(C), Ti{C). 
For later use, we need to compute some of the commutators of 0-© more in detail. 



-ih B ,e&d\ t . 
we also shall 



[At,!Rj£, = i(Ejj — Eu) — < 



For (j£) ? (pq): 



j -,,3-1 -r ^p= 3 
— A\),t-l 



'^A(,^_i ifj>land£>2 
if j = 1 and i > 2 
if 1 = 2 

(11) 



[At,K,j'£ 5 A{,3,. 

[At,a,j<, A(,3, : 



(— 5lp\t y SRJ q — SjpXe^Jq + SjqXt^,lp + $lqM,$t,jp) 



(12) 



B. Unforced equation 

For pure states in an orthonormal basis, the eigenvectors \ipj) — ej of the stationary Shrodinger equation 
are mapped into the diagonal density operator = Ejj- More generally, for quantum ensembles, 

after a suitable diagonalization, p — diag(iui, . . . , wn), ^ Wj ^ 1, J2j=i w j = 1- The eigenvalues Wj 
represent the populations of the various energy levels and provide a complete set of invariants for J5J, call it 
J = {wx, . . . , wn}, since © is isospectral. 

Proposition 1 Consider the system IjlOJI with Ha strongly regular. The state g is an equilibrium point of 
(llUf) for u — if and only if p — qqXq + q§ Af, . Furthermore, if gt ^ 0, then for u = 

1. u(p(0)) = h(p(t)); 



2 - ei M ,ji + et, %ji = const; 
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3. for St small, Qi,!R,je(t) ^ Qt,9t,jt(t + St) and Qt,%je(t) ^ Qt,<n,jt(t + St) V (jt) £ Ft{p). 
Proof. When u = 0, for a given p = g \ Q + p^ + pt, 

-i[H A , pt,} = (13) 
-i[H A> p] = -i[H A , pt] = -i[H A , Qt,m,je^iMje + Qt,3,ji^t,%je] 

— (£j ~ £l) {QiM,jiM,?s,ji ~ Qt&jiM&jt) ■ (14) 

In terms of the coherence vector g and using the isomorphism (|A1(I (meaning Ag(t) ~ — i[H A , p]), from 
(|13|l if Qt = 0, Ag = 0, i.e., p — qq\q + pt, is a fixed point. To show the other direction, notice that 
in (|14fl £j — £( ^ V (j£) l^j<£^N — 1, since H A is nondegenerate. Hence whenever Tt(p) ^ 0, 
Aq = Aqi ^ 0, because of the invariance of the t pq subspaces under f), see also (0). Therefore when Qt ^= 
the unforced system flows along nontrivial periodic orbits. Condition^of the last part also follows from 0. 
Since Q^(t) = £>ij(0), it must be ||f?e(i)|| = const Vt. This, together with the invariance property Q) yields 
Finally, Item [3] follows from pt never being fixed under the flow of the drift. To see it, consider a small 
time increment St. In the first order approximation, one can write 

g(t + St) = e StA Q {t) = (I + StA) g(t) 

i.e., the increment at St is given by l|14|l and the claim follows from the fact that H A is strongly regular, i.e., 
£ 3 -£l^£ p -£ q ,{jl)^{pq) J^i.p^q. " □ 

Since dim(j4) = n > N , the stationary Liouville equation has more eigenvalues than those referable to 
the eigenvalues of the corresponding Schrodinger equation. From A — —ia,dn A , these are the roots of the 
Lie algebra su(N) computed at the element H A of the Cartan subalgebra f), and, from i|14fl . they correspond 
to the transition frequencies £j — £i. The regularity of H A guarantees that these extra eigenvalues of A 
are all nonzero: dim (ker(j4)) = N — 1 = dim(t)), thus providing an alternative proof of the first part of 
Proposition 

For tr(p 2 ) = 1, the mapping — > |^)('(/'| eliminates the ambiguity in the (unobservable) global phase: 
= e lv \ip)(ip\e~ lv y tp £ R and, as before, the same property holds also for mixed states. Proposition^ 
affirms that, consequently, the corresponding one-parameter orbit passing through each eigenstate (due 
to the global phase) collapses into a fixed point of the unforced Liouville equation. Rephrasing in terms of 
density operators (part of) Proposition^ we have the following. 

Corollary 1 Any diagonal density operator is a fixed point of J5J when u — 0. More generally, for any 
density operator both the diagonal part and the trace square norm of the off-diagonal part are integrals of 
motion of J5J when u = 0. 



III. STRUCTURE OF THE STATE SPACE: COMPLEX FLAG MANIFOLDS 

It is possible to give a more thorough interpretation of Proposition ^ by studying the structure of the 
manifold in which p is living, call it S. S is a connected, simply connected submanifold of S"" 1 (the {n— 1)- 
dimcnsional sphere of radius ||f?||) whose dimension depends on the multiplicities of the eigenvalues of p. For 
N > 2, S C S n_1 , since the Lie group exp (ad^N)) is not acting transitively on the entire § n_1 . S, instead, 
is a homogeneous space of exp(ad 5u (Ar)), the action being left matrix multiplication, and can be described 
as a (co)adjoint orbit of SU(N) on its Lie algebra as follows. Consider a diagonal density Qt, £ S. Call 
Cg h the stabilizer of g^, Cg h = {g £ exp (ad 5u (jv)) s.t. gg^ — g^}. Because of the identification (up to the 
imaginary unit) of the density operators with (a convex set in) su(iV), the coset space exp (ad 5u (7v)) /Cg. 
is the adjoint orbit of SU(N) on its Lie algebra passing through g^. Because of transitivity, this orbit 
can be identified with S: S — exp (ad au (jv)j /Cg h Qij- The dim(5) is always even (each (co)adjoint orbit 
has a symplectic structure as is well-known). The orbit S is transverse to f) and meets f) in a number of 
disjoint points equal to the number of distinct permutations of the entries of g^. Such number is equal to 
the cardinality of the Weyl group as well as to the Euler characteristic x(<S) of the orbit, see [rH l39l | and 



G 



Theorem E.2 of [T^. These points form the vertices of a polygon in the N — 1-dimensional eigenensemble 
sitting in () and are sometimes denoted Weyl chambers. Inspired by the § 2 case (see Example below) , 
we shall call them antipodal. If = diag (wi, . . . , Wn), Ylj—i Wj = 1, ^ Wj 1, then the x(<5>) — 1 
antipodal points are given by diag . . . , m^at)) with er(l), . . . , o~(N) a permutation of 1, . ..,N such 

that diag («) ff m, . . . , Wo-(n)) 7^ Pf;- While the topology of the diagonal coset representatives is particularly 
easy to visualize, the entire orbit enjoys the same topological structure of g^. To see it, simply notice that 
applying a rotation in exp(ad 5u ( A r- ) ) to two or more diagonal antipodal points they remain antipodal. Since 
exp(ad su (jv)) ac ts transitively on 5, this is true on the entire orbit: each g £ S has x(^) — 1 antipodal 
states in S. What is not known a priori is the isotropy subgroup, which depends on J . In fact, the convex 
set of iV-level density operators is foliated into leaves of different dimensions, depending on the number of 
distinct Wj and on their multiplicities. For example, for a pure/pseudopure state J = {w±,W2, ■ ■ ■ ,1^2}, 
wi 7^ w 2 , Wi + (N — l)w 2 = 1, S = exp(ad 5U ( A r))/ (S 1 x exp(ad su ( A r-i))) and dim(<S) = 2N — 2. In the 
pure state case J — {1,0, ...,0}, it is well-known that the map \ip) — > \i/^)(ip\ can be seen as a Hopf 

fibration: § 2W_1 > S = CP 1 ^^ 1 , with fibers representing the global phase. At the other extreme, if 
J = {wi,W2, ■ ■ ■ , wn}, Wj ^ we, Ylj=i w j = 1; then S — exp(ad su (jv))/ (S 1 ) of dimension TV 2 — N. Hence 
if to = dim(«S), 2N — 2 ^ m ^ iV 2 — N, to even. In between lies the flag manifolds with flag determined by 
the multiplicities of Wj. If such multiplicities are given by ji, . . . , jj>, j\ + . . . + jp = N, 2 ^ £ ^ N, 

S = exp(ad su(A r))/ (exp(ad 5u(il) ) x . . . x exp(ad su(if) ) x (S 1 )^ ^ . 

When £ = 2 we have Grassmannian manifolds. Normally, in the literature these are known as complex flag 
manifolds and are given directly in terms of unitary group actions, see [ll llll 39] : 

S = U(N)/ (U(ji) x ... x U(u)) , jx + ■■■+ jt = N, 2^£^N. 
= SU(N)/ (sufa) x ... x SU(j t ) x (S 1 )'" 1 ) 

In terms of unitary actions, the two extreme cases of pure states and all different eigenvalues are, respectively, 
S = U(N)/(U(N - 1) x £/(l)) = SU(N)/ (SU(N - 1) x S 1 ) and S = U(N)/ (U(1)) N = SU(N)/ (S 1 )^ -1 . 
The description adopted here is just an isomorphic real representation of such complex flag manifolds deriving 
from the use of the adjoint representation. 

Example 1 N = 2, J = {1,0}. The case N = 2 is the only easy one, as S = S 2 ~ CP 1 . On the great 
horizontal circle of S 2 , g = gt- In terms of the Bloch vector, the diagonal antipodal states become the north 
and south poles of the Bloch sphere, 



p Vl = |0)(0| =diag(l, 0) ^ g Vl 
Pv 2 = |1}(1| = diag(0, 1) g V2 



1 l T 
° 75 







1 

V2 



T 



and f), dim(f)) = 1, corresponds to the vertical line passing through g Vl , g V2 . Everything extends unchanged 
to mixed states, since S is still equal to S 2 regardless of the purity. Since each S crosses f) exactly twice, 
x(S) — 2. For any g £ S the antipodal state is — g. 

□ 



Example 2 N = 3, J = {1,0,0}. Since the isotropy subgroup in this case is 5*0(3) x S 1 of dimension 4 
(recall that dim(ad su (3)) = 8), dim(6>) = 4 and x(«S) = 3. Following the standard ordering convention, the 
3-level Gell-Mann basis (see e.g. |2j|, p. 99) is 

{At,SR,12, A{ i Q i i2, A^i, A{ ; sjf,13j A{ : Q : i3, \i,SR,23> M,$l,23i ^^2}. 
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The three diagonal antipodal states are 

p Vl = diag(l, 0, 0) g Vl = 

Pv 2 = diag(0, 1, 0) <^ q V2 = 

Pv 3 = diag(0, 0, 1) q V3 = 



71 To 
"A 776 
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The structure of S C S 7 is studied in detail in 0, El EH HI] ■ In terms of the coherence vector, one has that 
each state of a triplet of antipodal states is at an angle of from the other two, see Fig. In particular, 
if Q G S, then —q^S: the single antipodal point of the case N = 2 is replaced by two symmetrically 
distributed and equidistant antipodal points. As expected, only A^i and Af,^ are of concern when p is 
diagonal. In an attempt to visualize the entire S manifold, one should replace each edge connecting to 
vertices with a sphere S 2 . The same numbers occur for pseudopure states, which have the same twofold 
degeneracy and are obtained by rescaling down the coherence vector g Vj by a constant factor. For the all 
different eigenvalue case ptf — diag(iui, 11)2, W3), Wj ^ Wk, which is the generic case, the stabilizer is the 
torus S 1 x S 1 , S = exp(ad su (j V ))/(§ 1 x S 1 ) and dim(6>) = 6. The only diagonal matrices that are conjugate 
with p^ are its five element permutations, i.e., x(<S) = 6 in this case. The six vertices are given by 



diag (a, b, c) 
diag (6, a, c) 
diag (c, a, b) 
diag (c, b, a) 
diag (6, c, a) 
diag (a, c, b) 



> Qga 



0^ a+b-2c 

00^ a+ ^ 2c 

00^ a+ ^ 2b 

00^ b+ ^ 2a 

00^ b+ ^ 2a 

0^ 2d^» 



(15) 



The 6 vertices correspond to -y rotations, plus their reflections around the three axes bisecting the triangle 
of pure states, see Fig. Represented in the (Af,,i, Af,^) plane, the eigenensemble of p has the shape of a 
hexagon inscribed in the triangle having as vertices the pure states, see Fig. ^ 
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FIG. 1: Vertices of the eigenensemble in the (A^i, plane for N = 3. Vertices joined by the solid line: pure 

states. Vertices joined by the dotted line: "generic" mixed state with all different eigenvalues: a — 3/5, b = 3/10 
and c = 1/10. Vertices joined by the dashed-dotted line: another "generic" mixed state with eigenvalues a — 9/10, 
b = 2/30 and c = 1/30. 
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IV. FEEDBACK STABILIZATION FOR iV-LEVEL QUANTUM ENSEMBLES 

A. Problem formulation 

For the system 111 Oil , we are interested in the problem of tracking a periodic orbit. More precisely, the 
stabilization problem is the following. 

Given p d € S, find u — u(g d , g) such that, for t — > oo, g — > g d , where 

g d = A d g d , A d = h Ad ■ ad A „ • (16) 

This is a full state tracking problem which, from Proposition ^ reduces to stabilization to an equilibrium 
point when p d = g X + g dA ■ Af,. 

B. A modified Jurjevic-Quinn condition and antipodal points 



for h/>) ( 
24]. It 

v = U 

V = V(g d , g) = -((A d g d , g}} - ((g d , Ag)) - u((g d , Bg 



The algorithm for the feedback design resembles the one used for h/>) discussed in |13,|21|,|28|,|29, 36] and 
indeed the standard Jurdjevic-Quinn method for bilinear systems [24|. It consists in choosing a distance-like 
candidate Lyapunov function V — V(g d , g). From JJjJ, consider 1 V — \\g\\ 2 — {(g d , £»))■ If g d obeys to JT^J, 



= ({g d , (A d - A)g)) - u{(g d , Bg)). (17) 

This expression has a drift term which is in general sign indefinite and for which there is no global smooth 
feedback compensation. Whenever V can be rendered homogeneous in u (namely when Ha a = Ha, which 
is assumed thereafter), there is an obvious choice of feedback guaranteeing positive semidefiniteness. For 
example, 

u=((g d ,Bg)) (18) 

makes V = —v? ^ 0. In order to apply LaSalle invariance principle, one can try to adapt the "ad-condition" 
of Jurdjevic-Quinn |24( to the case at hand. The largest invariant set £ in J\f = {g s.t. V — 0} can be 
computed imposing u = % = . . . = = 0, where, in correspondence of u = 0: 



and, similarly, 



d ^ = {~lYe T d [A,...,[A,B]...]g = Q. (19) 



For \ip) which is an eigenfunction, this condition, used implicitly in 36], is made explicit in |2£| as a Kalman 
rank condition on the linearized tangent system. It can be reformulated for densities as follows. For a given 
g d G S, write g = g d + Sg. The linearization of l(TtJ|) around g d 

^- = A5g + bu (20) 
at 



1 Notice that rather than the "distance inherited from the S" — 1 sphere" V used here, one could use the 1" distance ||(5g|| 2 = 
llffd — q\\ 2 33 candidate Lyapunov function. Up to a scalar factor, the two give the same gradient, hence the same control 
design. 
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where 

b = Bg d , (21) 

is living on T Qd S. Since dim(T etJ 5) = m, if the Kalman rank condition 

rank [b Ab . . . A m_1 6] = m (22) 

is satisfied, then in the same spirit of the original Jurdjevic-Quinn work, this implies that £ contains no 
other trajectory than Qd, at least locally. Of course if pd = Qq\q + Qd,t) ■ A, then under (|16|l p d (t) = Pd(0) and 
6(t) = 6(0) in J2U- However, when pd(0) is not diagonal, Pd{t) is time- varying and so is b(t), complicating 
the verification of (|22|l considerably. Furthermore, in this case the condition ill'l't does not have a global 
character because of the topological structure of S. To see this, consider the case of diagonal density operators 
in which the linearization (|20|l is time-invariant. Call g p an antipodal state of Qd- Then also g p is diagonal 
and so is 6g p — g p — Qd- Hence A8g p = 0, i.e., (121 )[) has vanishing drift in correspondence of the x(S) ~ 1 
antipodal states of gd- Therefore, in spite of the Kalman condition (|22(l . in this case when checking LaSallc 

invariance principle = for H20|l since u = in A/". This argument can be generalized as follows. 



Proposition 2 Given gd € 5, any o/ i/ie other x(<S) — 1 antipodal states g p € S is an equilibrium point for 
the system (|10J) urai/i £/ie feedback Ijl8(l . 

Proof. We make use of the isomorphism (|A1I) . For any given diagonal, [— iHB,p P ] £ * s 

diagonal i.e., Bg p = g' p V But if p d diagonal p d = g Q X + g t)d ■ X then u = ((g d , Bg p )) = ((g d , g' p ^)) = 0, 
since f) _L t. Hence no feedback is produced. When instead g p , gd & S are antipodal but not diagonal, 
then (by construction) 3 R G exp(ad 5U (7v)) such that g p = Rg p and g~d — Rgd are both diagonal. The 
skew symmetric matrix RF (—i&dH B )R belongs to the adjoint orbit of exp(ad 5U (jv)) in ad su (jv)j hence 3 
C £ su{N) such that R T BR = ad c . Therefore u = ((g d , Bg p )) = ((g d , R T BRg p )) = ((g d ,&d c g p )} = 
because adc£> P — [C, p p ] £ t/i, while p d £ i). □ 

Corollary 2 For pure or pseudopure states, the N antipodal points of S are all equidistant. For pure states 
they are also maximally distant. 

Proof. It is enough to notice that for any triple of antipodal points p Pl , p P2 and p Pa (p d included), 
tr {p Pl p P2 ) = tr (p Pl p Pa ) = tr (p P2 p P3 ), hence V(g Pj , g pt ) = const V j, £ = 1, . . . , N, j ^ I. For pure states, in 
addition, tr (p Pj p Pe ) = 0, hence V(g Pj , g Pe ) = ti(p pj ) are maximally distant. □ 

As will be shown in next Section, the antipodal points are not the only states lacking attractivity, and the 
linearization alone is not enough to investigate the domain of attraction of the feedback stabilizer. 



Remark 1 The trajectory tracking problem presented above admits a reformulation as a point stabilization 
for a nonautonomous system. Consider a frame rotating with A. Call g d and g the new reference and state 
vectors. Then g(t) — e~ tA g(t) and gd{t) — e~ tA gd{t) — ffd(0), i.e., the reference trajectory becomes a fixed 
point. Using a variation of constants formula, we obtain for i|l(J|) 

= ue- tA Be tA g 

= e(o)- { } 

The Lyapunov distance is V = \\dd\\ — ((dd,gj) and its derivative V = u((gd,gj). The uniformity of the 
asymptotic stability for the nonautonomous system l|23|l with the same feedback stabilizer as l|18|) follows 
directly. 
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C. Time-independent convergence conditions 

To formulate a convergence condition in a more geometric manner, rewrite (|19fl in terms of bilinear forms 
of skew-symmetric matrices as follows. Call 

2XT = span J B, [A, B] , . . . , [A,...,[A, B] .] I 

v a times / 

and 20^4 = spanjA, 2H Q }. W a contains skew-symmetric matrices, and the conditions ^ — 0, . . . , a, 
written compactly as gJW a g, are bilinear forms. If q// Qd then g^W a g = 0. On a sphere q// Qd means 
g = ±Qd- However S is only a submanifold of S" _1 , and —Qd may not belong to it at all. 2U" is invariant 
to the so-called "ad-brackets" but not necessarily a Lie subalgebra. Of course if for the system p0(l the 
ad-brackets are generating, i.e., if, for some a, = ad 5u (jy), then almost global convergence on S is always 
verified. However we have the following negative result. 

Lemma 1 If N ^ 3, for the system (|10J) with Ha strongly regular and Graph(_ffs) connected, Lie(2H^) = 
ad 5u(A r) but W\ C ad su(Af ) Va>0. 

Proof. Since Ha is strongly regular and Graph(i/s) is connected, it follows from Theorem 2 of that the 
smallest subalgebra containing —iHa-, —iHb is s\x(N). Hence the same holds for the adjoint representation. 
For the second part, recall that dim(f)) — N — 1. From the Lie bracket relations J7|), A E adf,, B 6 ad{ 
implies [A, . . . , [A, B] . . .] e ad{. Even adding A, W\ alone cannot fully generate adf, for any a. □ 

The first part of Lemma^is also known as the strong accessibility condition |30j . Since ad 5u (jv) is compact, 
it suffices for controllability. The second part is the Jurdjevic-Quinn condition mentioned above. If Hb is 
not off-diagonal, then the statement of LemmaQ] should be reformulated as "N > 3". 

The following Theorem provides a time-independent condition for asymptotic stabilizability to any p c i G S, 
and a global description of the region of attraction of the controller. 

Theorem 1 Consider the system <|10|) with the feedback (|18fl . where pd G S obeys to (|16f) . Assume that Ha 
is strongly regular and that Hb is such that (/i_b,k,j <j+Xi 7^ (0, 0). An initial condition p(0) G S is 

asymptotically converging to Pd{t) if 

1. p(0) is not an antipodal point of pd(0), 

2. T([H B , p d ])nT(p(0))^0, 

3. CardJ- t ([F S) Pd ]) ^ m /2 

where Card denotes the number of pairs of indexes in Ti . 

In order to prove the Theorem we need a few preliminary results. 

Lemma 2 Under the assumption of strong regularity of Ha, the following three conditions are equivalent: 

1. the Kalman rank condition (|22|) is satisfied; 

2. rank {W 1 ' 1 Q d ) = m; 

3. C&rdF t ([H B , P d}) >m/2. 

Proof. Given C £ su(7V), strong regularity of Ha implies that C, [Ha, C], . . . , [Ha, • ■ ■ , [Ha, C] . . .] are all 
linearly independent up to a number a — 1 , a — 2 CarcLF{(C), of nested Ha commutators, see Theorem 2 in 
0- Using C = c ■ X and the isomorphism given by the adjoint representation, the vectors c, Ac, . . . A a ~ 1 c 
are all linearly independent. If c = b = Bgd as in l|21|l . then this is the Kalman controllability condition 
provided a> m. If b = &(, + &{, A\ = 0, hence only the off-diagonal part of [Hb, Pd] matters. The support 
U{[Hb, Pd}) intersects a number of "root spaces" tji (each has real dimension 2) equal to Cax&J-t{[H b , Pd}) ■ 
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Furthermore, since —iH A G su(iV), the invariance property Q applies. Written in terms of the original 
commutators 

T t {[H B ,Pd\) = M[Ha, [H B ,p d ]}) =■■■ = T t {[H A , [H A , [H B ,p d ]} . . .]). (24) 
For the £-th order commutator of 25J a_1 , one has the binomial-like expansion: 

£ times 

= A t B + {-l) 1 (^jA l - 1 BA + ... (25) 

+ (^[^ABA 1 - 1 + {-ifBA 1 . 

After the linearization around pd, only the first term of this expression is retained. Since —iH A £ f), from 
0, f { (A e g d ) = ft (fid), while ft, = 0, VI ^ 1. From (J23J, .F{ is the same for all terms in and 

similarly, 

^, (Bf> d ) = ([A, = . . . = ? t ([A, . . . , [A, B] . . ]g d ) . (26) 

In summary, strong regularity of H A guarantees the full spanning of a linear space whose dimension is 
determined uniquely by CardJFf ( [H B , Pd])- This space is identifiable with the tangent space T Qd S as well as 
with W m ~ 1 Qd- The equivalence of the three conditions follows consequently. □ 



Remark 2 Lemma^and Condition [2] of Lemma[2]imply that although the vector space 2U" l_1 is never the 
entire Lie algebra ad su (jv) acting transitively on <S, it may nevertheless span the entire tangent space at a 
point. The same holds for the Kalman controllability. 

Remark 3 In general CaidT^H B ) ^ Cax&Tt^H B , Pd]), hence the controllability of the linearization de- 
pends from the reference trajectory pd chosen. The meaning of Lemma [5] is that in order to have linear 
controllability the reference trajectory (|16J) must be "rich enough" along H B . 

Remark 4 While the Kalman condition i122[| seems time-varying as soon as pd.t ^ 0, the equivalent condi- 
tion of Lemma [21 is always time-independent since J-([H B , pd(0)]) = J-([H B , Pd(t)]). 

Remark 5 The conditions of Lemma [3 depend on pd, H A and H B but not on the state p, meaning that 
alone they are not enough to guarantee convergence of a given p(0). 

The Lyapunov derivative in l|17fl is made homogeneous in u by the cancellation of the drift term and 
therefore the notion of attractivity provided by V must be rendered invariant under such flow (in a way 
similarly to the orbital stabilization problem, see Q). The following Lemma gives an alternative attractivity 
condition which is fully invariant under the drift and generically (i.e., almost always under e tA ) equivalent 
to the usual Lyapunov convergence property. This last in fact may fail in isolated points: certain critical 
points of V are not invariant under the flow of the drift (see Section Ivf . 

Lemma 3 Consider the system (|1U|) with the feedback 118fl . where pd obeys to (|16fl . If H A is strongly regular, 
the following conditions are generically equivalent under the flow of the drift term: 

1. T([H Bl p d })nT(p)^0; 



2- V(g d , q) < 0; 

3. g TW a Q = {z , Zl ,...z a }, Zj ^0. 
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Proof. Clearly V = — ((gd,Bg)) 2 < implies T ([Hb, Pd]) H.F(/o) 7^ 0. To prove that also the contrary 
is generically true, it is enough to show that when T ([Hb, Pd]) H T (p) 7^ the zero crossing of the inner 
product can occur only at isolated points along the trajectories of the closed loop system. Assume Item^ 
holds and, at time t, ((g d ,Bg)) = 0. If St is a small time increment, then from Item ^ of Proposition ^ 
J- ([Hb, Pd}) and T (p) remain the same, while, from Item|31of Proposition^ 

(Qt,$t,jt(t + St), g^je(t + St)) ^= (ge,x,ji(t), gt,Q,ji(t)) 
Qf>,j(t + 8t) = Q^j(t) 
(Qd,t,n,jt(t + St), Qd,t,^,ji(t + St)) 7^ (Qd,t,»jt(t), Qd,t,S,jt(t)) 
Qd,f),j(t + St) = Qd,t),j(t)- 

If Ttj ([Hb, Pd]) H Tt, (p) 7^ 0, then from the last row of (JSJ only g^t matters in the computation of 
Tf, ([Hb, Pd]), and ((gd,t(t + St),Bg(,(t + St))) ^ since g djt (t + St) ^ Qd,t(t), while g h (t + St) = g h (t). 
If, instead, Tt ([Hb, Pd]) H Tt (p) 7^ 0, then we have two possible contributions to consider: Tt ([Hb, Pd,t)]) 
and Tt ([Hb, Pd.e])- In the first case the conclusion follows from the same argument used above since now 
Qdfy(t + St) = g d fy(t) while gt(t + St) 7^ gt(t)- In the second case it follows from the observation that 
Tt ([Hb, Pd,i]) H Tt (p) 7^ implies Tt(pd,t) 7^ Tt (p) (see the explicit computations of the commutators in 
(23 )■ The general case T ([Hb, Pd])^T (p) 7^ is the sum of the two situations just described. Concerning 
Item |3 it is enough to notice that generically ((gd, Bg)) 7^ if and only if ((gd, [A, . . . , [A, B] . . .]g)) 7^ 0. 
The argument is of the same type used in the proof of Lemma|21 For example if Tt ([Hb, Pd,t\) H Tt (p) 7^ 
then just apply If, instead, we are in the case T$ ([Hb, Pd,t\) H Tt, (p) 7^ 0, then the only useful term 

in the expansion (|25|l is the last one, but this is enough to prove the claim. The genericity of the argument 
can be shown as above. □ 

Proof, (of Theorem^ Consider the set N '. We want to show that the largest invariant set £ in TV is given by 
pd only. Condition [3 guarantees that locally around pd(t) there is no other equilibrium point in Af, as, from 
Lemma |2 the linearization at pd is controllable. Hence pd is a locally asymptotically stable equilibrium for 
the closed loop system and pd is isolated in Af. Consider p e £ Af, p e 7^ pd- This implies p e disjoint from p d and 
y(Qd, £>e) > 0. We need to show that p e must be a repulsive equilibrium for the closed loop system 2 . For g e 
which is an antipodal point of a pure state g d , this is follows from Corollary|21 since V(f?d> Qe) is maximal in 
S while V ^ 0. For any other p e 6 Af, it is enough to perturb p e to p e £ S so that T ([Hb, pd]) H T (p e ) 7^ 0. 
It is always possible to do this in a neighborhood of p e since T([Hb, Pd]) has cardinality at least to/2 
and (h,B,gi,j j+i, h>B,$s,j j+i) 7^ (0, 0) implies that Graph(iJe) is connected and that there is no subspace tji 
invariant under B. But then, from LemmaEl ((gd, Bg e )) 7^ and V(gd, g e ) < V(gd, g e ), i.e., p e is attracted 
to pd- To show that V(g e , g e ) increases, assume by contradiction that 

V(g e , Qe) = ((Qe, Q e )) + ({Qe, he)) = -{{Od, B Q e )) «f?e, B Q e )) < 0. (27) 

Consider the geodesic line in S connecting gd with g e : g<f,(s) = gd + <fi(s) such that cp(0) — and cp(s e ) — 
ge — gd- Along this line, 

V(gt(s),g e ) = -{(g d ,Bg e )) 2 - {(g d ,Bg e ))({(l>(s),Bge)), s£ [0,s e ] 
is a function linear in (p(s) and such that, by the assumption l|27|) . 

V(Qt(0),Q e ) = V(g d ,ge) = -((g d ,Bg e )) 2 <0 
V(g4,(se),g e ) = V(g e ,g e ) < 0. 

But then V(g<p(s), g e ) < V s £ [0, s e ], V(g<f,(s), g$(s)) — 0, meaning that g e is attracted to the entire 
geodesic segment g t p(s), s £ [0, s e ], which is a contradiction, since g d is an isolated equilibrium point. Hence 
it must be V(g e , g e ) ^ i.e., g e is a repulsive equilibrium point. Therefore g e cannot belong to £ . From 
LemmaOll all conditions i|19|) are satisfied or violated simultaneously respectively when V = or V < 0, i.e., 



2 Since Q e may not be isolated in Af, the term repulsive has to be intended as "semi-repulsive" 
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when !F([Hb, Pd]) H F {p e ) = or ^ 0. Hence outside M the Jurdjevic-Quinn condition applies and p(0) 
must converge to Pd(t) since any other p e G TV is repulsive. □ 

Remark 6 Condition|21of Thcorcm^s obviously a necessary condition for convergence. Condition|3instead 
is sufficient but not necessary, see Example 1 in Section Ivl 

While, from Lemma [21 the linear span at Qd of the linearized system and of the W a yield a space of 
the same dimension, Item of Lemma 01 holds for the bilinear forms but it is in general not true for the 
linearization. 

Corollary 3 For Ha strongly regular: 

1. rankSm Q - 1 £ d = rank[bA6...A"- 1 6] , V a = 0, . . . , m - I; 

2. ((g d , i^JA, B] . . .]„)) ? <^((g, A a b)) ± 

a times 

Proof. The first point follows from the strong regularity of Ha and from (|26() . which implies that the 
maximum number of independent vectors in the two sequences above is the same for all a. The second from 
(25) and A a gt, =0. □ 

The consequence is that the linearization alone is inconclusive about the region of attraction of the reference 
trajectory in the closed loop system, while instead the ad-commutators completely specify it. 

Corollary 4 When Card^t ([-ffs, Pd]) TO /2, the region of attraction of the system l|l(J|) with the feedback 
(JTSJ is given by1l = S\Af. 



D. Global stabilization and topological obstructions 

The notions from differential topology used in this Section are recalled in Appendix^] A compact manifold 
like S cannot be globally asymptotically stabilized because it lacks the contractivity property, i.e., it is not 
homotopy equivalent to a point, see 0, Proposition 1 and Theorem 1, and [33 . Proposition [21 suggests that 
for S this is due to the antipodal points. 

Proposition 3 For S , the x{S) ~ 1 antipodal points are irremovable topological obstruction to global stabi- 
lizability by smooth feedback. 

Proof. Contractivity is a necessary condition for global asymptotic stabilizability. For example, that S = § 2 
with a point removed is homeomorphic (and hence homotopy equivalent) to R 2 is well-known through the 
stereographic projection (see e.g. @, P- 34). Since it is known that the domain of attraction of an 
asymptotically stable point must be homotopy equivalent to R m for some m (3?|], then for TV = 2 this is 
enough to affirm that convergence can be rendered global up to the antipodal point. For N > 2, in order to 
show that the antipodal points are all obstructions to contractivity of 5, consider the equilibrium pd and one 
of its antipodal points p Pl . By suitable change of basis, pd and p Pl can be rendered diagonal simultaneously. 
By the trans versality of the coadjoint orbit on f), it is possible to determine a submanifold of S connecting 
Pd and p Pl and not passing through any other of the antipodal points. To see it, notice that for N > 2 it 
is always possible to adjust the basis ©-© so as to attain a 3-dimensional simple subalgebra of su(A^), 
which, as described in Section ITTll draws an § 2 orbit under the adjoint action (see [33] for the details of this 
construction). This is a well-defined compact submanifold of S and it is not contractible for what said above 
for N = 2. 

□ 
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V. A FEW CASES OF PHYSICAL INTEREST 



The methods developed above yield considerable insight into the stabilizability and convergence properties 
of a quantum density operator. A few interesting cases for TV-level systems are now described. It is followed 
by a more detailed description for systems with N = 2, 3. 

• Since m ^ TV 2 — N, and CardJ^ ([Hb,Pci]) ^ {N 2 — N)/2, (i.e., the maximal number of off-diagonal 
terms), each complex flag manifold S may admit a controllable linearization (depending on pd). 

• The assumption of direct coupling between nearest energy levels (h,B,m.,j i+i> hB,$s,jj+i) ^ (0, 0), is 
needed in order to exclude the existence of subset of S which remains invariant under the closed loop 
dynamics. It is a common assumption in most practical cases (dipole approximation |Tfijh See also 
Example 2 below (last item). 

• The full connectivity of Graph(Zfg) is neither a sufficient nor a necessary condition for asymptotic 
stability. 

• If pd is an eigenstate and p another eigenstate then there is never convergence, not even if Graph(_ffe) 
is fully connected, because p is antipodal to pd- 

• For pure states and not fully connected Graph(iJs), certain eigenstates are easier to stabilize than 
others. The easiest is the one of energy Sj such that the index j appears more often in Ti ([Hb, Pd])- 
In Example 2 below with Hb in J3SJ), it is easier to stabilize to the eigenstate of intermediate energy 
than to the ground state or to the most excited state. When Graph(iJe) is fully connected, there 
is no such difference. From Theorem ^ this does not mean that all initial conditions have the same 
convergence properties to a given pd- 

• If pd and p(0) are both block diagonal and the blocks do not overlap 



* . . . * 



* ... * 

then 



[H B , Pd] = 



* ... * 






* ... * 


* . 


. . * 


* . . . * 


. 


. 


* . . . * 


. 


• o. 



which implies J 7 ([Hb, Pd])H T (p(0)) — and V = 0, i.e., p(0) is not attracted to pa- Since tr (ppd) = 0, 
pd and p are as distant as antipodal states. 

• Not all states in Af are maximally distant from pd- Assume pd, p such that ([Hb, pd])^^ (/°(0)) =/= 0, 
Ft (Pd) — 0, J-( (Hb)^Fi (p) = 0. Also in this case T ([Hb, pd])^? (p(0)) = and p is not converging. 
However, since tr (pdP) 7^ 0, pd and p are not maximally distant. 

• A typical example of initial condition such that V(Qd(0), f?(0)) = but not invariant under the drift 
(see paragraph before Lemma OJ) is attained when Tt ([Hb, Pd]) H T\ (p(0)) ^ but [Hb, Pd(0)], p(0) 
both real or purely imaginary. This follows from Proposition ^ 



Pd 



* . . . * 



, P(0) 
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Example 1 (cont'd) Assume 



H A = 



V2 



1 

-1 



and H B = 



l B,tM,12 



V2 



1 

1 



Then 



A = -iv2/iA,6,iadA[, >1 = 2h A ,t,,i 



0-10 
1 




and B = -iV2h B ,t,3i,i2&d\ l3il2 = 2h B ,tM,w 




0-1 
1 



From Proposition ^ both ||f?{|| and are integrals of motion of the unforced dynamics, while the two 
components of Qt evolve according to a sinusoidal law. When applying Theorem ^ to the system plus the 
feedback l|18|) . we have the following for the closed loop system: 

• any pd has a single antipodal point which also is an equilibrium; 

• if pa diagonal, Tt ([H B ,pd]) = {(12)}, the linearization is controllable and any nondiagonal p satisfies 
Theorem ^ Hence any p(0) such that pt(0) ^ is attracted to pd diagonal; 

• if pd off-diagonal, CardJFj ([H B ,pd]) = 0, and the sufficient condition of Theorem ^ does not apply. 
However, ([H B ,Pd]) ^ and as long as ([H B ,Pd]) H Ft, (p(0)) ^ 0, i.e., whenever pt, ^ 0, p -> p d . 
This is a special situation due to dim(f)) = 1, and has no counterpart for N > 2. 

In summary, there is always almost global convergence except when Qd.^ = Qtj = 0, i.e., except when both 
Qd and g belong to great horizontal circles. □ 

Example 2 (cont'd) The drift of the system is given by 



H A = 



7T 



1 
0-10 




h A . 



(,.2 



V6 



1 
1 

0-2 



= hA,tj,i\,i + h A ,t),2^t),2- 



We shall consider the following control vector field 



Hr 



1 

71 



/lB,t,Sf,12 

/lS,{,5R,12 /l.B,t,St,23 

h B fW,23 



= h B {M,12Xt,?H,12 + h B .tM,23^tM,23, 



(28) 



which has Fe(H B ) = {(12), (23)} or, alternatively, 



1 

V2 



h>B,t,9t,12 h B ,t,9t, 13 

h B ,t,m,12 /lB,t,ffi,23 

^B,( : 3i : 13 h B ^Ma3 



B,{,SR,12At 12 



JB,e,5R,13At,5R,13 



iB,e,5R,23At,3f,23 



(29) 



which has a "fully connected" graph, F t {H B ) = {(12), (13), (23)}. 
A list of interesting cases is the following: 

• any of the (two for pure/psudopure, five for the generic case) antipodal points of any pd is also an 
equilibrium. 

• pd diagonal: only the off-diagonal part of p matters 

- p d pure (or pseudopure), e.g. p d = p Vl 

* H B given in (|28fl : Ft ([Hb, Pd]) = {(12)} =>■ the linearization is never controllable since 
2CardJF t (\H B ,p d ]) < 4 = m, hence Theorem ^ does not apply. Unlike the N — 2 case, now 
in general p(0) -/>p d ; 
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* Hb given in Ft ([H B ,Pd]) = {(12), (13)} =$> the linearization is controllable. Any p(0) 
such that Ft (p(0)) n {(12), (13)} ^ is converging. However, if one considers the pure state 



P(0) 




1 1 
1 1 



then Ft (p) = {(23)}, implying V^(Q) = u = ({gd{0), Bg(0))) = 0, i.e., the system is not 
converging to pa in spite of the Kalman controllability condition on the linearization. Notice 
how for this example rank (2U 3 £ d ) = rank (W 3 g(0)) = 4, while qJ20 3 q(0) = {0, 0, 0, 0}. 

— pd pure (or pseudopure), but pd — p V2 

* Hb either (|28l) or l|29l) : Ft ([HB,Pd]) = {(12), (23)} =>■ the linearization is always control- 
lable. Any p(0) such that Ft (p(0)) H {(12), (23)} ^ is converging. 

— pd with all different eigenvalues, e.g. pd = p gi 

* Hb in (|28|l : Ft {[Hb, Pd]) = {(12), (23)} the linearization is never controllable since now 
m = 6; 

* Hb in l(2"§l) : Ft ([Hb, Pd]) = {(12), (13), (23)} ==>■ the linearization is always controllable. 
Any p(0) such that Ft (p(0j) ^ is converging, any p(0) such that Ft (p(0j) = is antipodal. 

Pd — Qo^o off-diagonal 

— Hb in (|28|l and Ft (pd) ^= Ft (Hb) ==>■ linearization is never controllable, hence Theorem 2] does 
not apply and in general p(0) -/*Pdif)) 

— Hb in (|28|l and Ft (pd) (Hb) ==>■ GaidFt ([Hs,Pd]) is at least 2, implying that the lineariza- 
tion is controllable at least for pure/pseudopure states; 

— if Ft (pd) H Ft (Hb) 0, then also F^ (p(0)) matters for the convergence, see (lllfl ; 

— if Ft (pd) H Ft (Hb) = then convergence depends only on Ft (p(0)) (plus controllability), see 

If the control Hamiltonian is Hb = ft-t.sR.12At.sR.12 + ftt,SR,i3At,SR,i3, i.e., direct coupling between £ 2 and 
£3 is missing, then the sufficient condition of Theorem^does not apply. Assume for example 



Pd 





* * 
* * 



p(0) = 



* * 

* * 




Then Ft ([Hb, Pd}) = {(12), (13)} and Ft ([H B ,Pd]) n F t (p(0)) = {(12)}. However, p(0) -/>p d {t). 



□ 
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APPENDIX A: ON THE ADJOINT REPRESENTATION 

A representation of a Lie algebra g on a vector space X is a mapping O : q — > fl'(^) which is a Lie algebra 
homomorphism, i.e., a map which 

1. is linear Q(a\Ai + a 2 A 2 ) = &i@(Ai) + a 2 Q(A 2 ), V A\, A 2 £ q and V ct\, a 2 in the field of X; 

2. preserves the Lie bracket Q([A 1 , A 2 \) = [9(Ai), <d(A 2 )}, V Ai, A 2 € g. 
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So a representation 9 assigns to each A £ g a linear operator 9(A) e Q^{X). A particularly useful 
representation is the adjoint representation on X = R". If for an n-dimensional Lie algebra g we choose the 
basis A\, . . . A n then the Lie brackets of the basis elements are [Aj, A k ] — c jk^- The components of 

the 3-tensor Cj k are called structure constants of the Lie algebra with respect to the basis Ax,.. . A n . The 
adjoint representation of g, ad B , with respect to the basis Ax,.. ■ A n is the representation having as basis 
elements the n x n matrices of structure constants Ax, ■ ■ ■ , A n , Aj = ad^. = [Aj, ■ ] of entries (Aj) gk — Cj k . 
Notice how the two free indexes k and £ identify respectively the columns and the rows of the new basis 
elements. 

In general, the adjoint representation of a linear Lie algebra is a derivation of the algebra and corresponds 
to the infinitesimal representation of all the one-parameter groups of automorphisms. For a semisimple 
compact Lie algebra g, the main features of ad B are (see e.g. (23, p. 39 and 129): 

• it is a real semisimple Lie algebra; 

• it is isomorphic to g; 

• VA, Beg: [ad^, ad s ] = ad[ A ,B]- 

Let us spend some more words on emphasizing how the "linearity" of the adjoint representation may be 
intended, which is one of the leitmotifs of the paper. If B € g has the expression B = bxAx + . . . b n A n , then 

r -1 T 

as long as we keep the basis fixed, B is uniquely identified by its vector of components: B ~ b = bx ■ ■ ■ b n 
Then 

n n 

[A, , B] = [A, , J2 h M - E $ ( &d ^) ek h = Ajb. (Al) 
fe=i fe, i=i 

We will often make the double substitution {^4, B} ~ {A, b} which will correspond to replacing the (bilinear) 
matrix commutator [ • , ■ ] : g x g — ► g with the linear operation ad B x R ra — » IR n , i.e., left matrix multiplication. 



APPENDIX B: A FEW FACTS FROM TOPOLOGY 



The material in this Appendix is taken from standard texts on (differential) topology e.g. [UH^]. Let X , y 
be topological spaces and f,g : X — > y be continuous maps. The mapping / is said homotopic to g if "it can 
be continuously deformed to g" , i.e., if 3 a continuous mapping h : X x [0, 1] — > y such that h(x, 0) = f(x), 
h(x, 1) = g(x) V x S X. X, y are said homotopy equivalent if 3 maps / : X — > y and g : y — > X such that 
g o f and fog are homotopic to the identity maps in X and y respectively, y is said contractible if the 
identity map on y is homotopic to the constant map y — * x a for any x £ y. A space is contractible if and 
only if it is homotopy equivalent to a point. No compact manifold is contractible. Homotopy equivalence 
is an equivalence relation on topological spaces and the classes of homotopy equivalent spaces are called 
homotopy types. Examples used in this paper are: 

• § m — {p} is homotopy equivalent to K m (p is any point of § m ) ; 

• R m — {0} is homotopy equivalent to S m_1 

• R m is homotopy equivalent to a point (and hence to any contractible space). 
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